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Abstract 

The Planck constant (h) plays a pivotal role in quantum physics. Historically, it has been 
proposed as postulate, part of a genius empirical relationship E = hu in order to explain the 
intensity spectrum of the blackbody radiation for which classical electrodynamic theory led to an 
unacceptable prediction: The ultraviolet catastrophe. While the usefulness of the Planck constant 
in various fields of physics is undisputed, its derivation (or lack of) remains unsatisfactory from 
a fundamental point of view. In this paper, the analysis of the blackbody problem is performed 
with a series expansion of the electromagnetic field in terms of TE, TM modes in a metallic 
cavity with small losses, that leads to developing the electromagnetic fields in a complete set of 
orthonormal functions. This expansion, based on coupled power theory, maintains both space 
and time together enabling modeling of the blackbody’s evolution toward equilibrium. Reaching 
equilibrium with a multimodal waveguide analysis brings into consideration the coupling between 
modes in addition to absorption and emission of radiation. The properties of the modes, such 
as spectral broadening, losses and lifetime, then progressively become independent of frequency 
and explains how equilibrium is allowed in good conductor metallic cavities. Based on the free 
electron relaxation time in gold, a value of h = 1.02 x 10 -34 J-s for the reduced Planck constant is 
found and the uncertainty principle is also emerging from this a priori classical study. The Planck 
constant is then obtained no longer as an ad hoc addition but as a natural consequence of the 
analysis taking boundary conditions into account as into optical resonators. That analysis based 
on finite-spacetime paradigm, also shine new light on the notion of decoherence in classical optics 
and electrodynamics. 
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I. INTRODUCTION 


The ninetieth century saw great progress in both electromagnetism and thermodynamics 
with radiative heat transfer studies positioned between both. The latter studies of light- 
matter interactions led Gustav Kirchhoff to introduce the concept of perfect black bodies and 
the universality of their spectral radiance as a function of wavelength and temperature [Tj. 
As light was identified with electromagnetic (EM) waves, the blackbody became an ideal 
tool to connect energy notions with the thermodynamics of electromagnetic radiation. The 
initial EM analysis of thermal blackbody radiation by Lord Rayleigh and J. Jeans is well 
known to have supported the introduction of a light quantum due to a discrepancy in the 
radiative energy density distribution at high frequencies, dubbed ultraviolet catastrophe. 
Great successes of the radiation quanta in explaining both the photoelectric effect and 
blackbody radiation led the EM field theory and the Maxwell’s equations to fade away from 
the spotlight. Over the years, semi-classical models without light quantization turned out 
to suffice in many cases of light-matter interaction, including the photoelectric effect [2) 
[3J. Besides, time-independent Schrodinger and Maxwell wave equations can both yield 
Helmholtz equations or other mathematically similar forms with a correspondence between 
potential and index of refraction. The reader can refer to [I] for a review on the wave- 
particle notions. Let us present a citation from the author (p.9), about the Klein Gordon 
equation of the relativistic quantum mechanics; ”An equation for a zero-mass particle of 
spin one (photon) in relativistic quantum mechanics turn out to be none other than the 
classical wave equation for the electromagnetic field of the 19 th century that predates both 
relativity and quantum physics”. Being inherently relativistic and adaptable to Hamiltonian 
formulation, the EM field theory was revisited with quantum electrodynamics (QED), a 
theory which explains with excellent precision experiments, as example, Lamb displacement, 
hydrogen hyperfine structure, and serves to analyze with success many others interactions. 
QED still however calls upon infinite plane wave approximations. An infinite plane wave 
basis used for the spatial dependence of the EM field is in contradiction with the finite 
size of a blackbody cavity, which incidentally should be seen as an optical resonator with 
boundary conditions allowing a penetration depth of the EM field in the walls. According 
to classical physics, the blackbody problem is analyzed under the hypothesis that atoms 
located on the cavity wall can only absorb and/or emit electromagnetic energy thus giving 
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rise at thermodynamic equilibrium to a system of stationnary waves (stationnary modes) 
in the cavity. Studying perfectly reflecting walls and a system at equilibrium lead to a 
reduction of the electomagnetic modal structure from E(r,t ) to g(t)E(r). It is thus not 
surprising that classical electromagnetism fails to explain how equilibrium is reached since it 
defines equilibrium as a time and a space separated functions and neglects a very important 
phenomenon of multimodal propagation in a physical medium, namely: mode coupling. 
What really leads to the ultraviolet catastrophe in the analysis of this problem is thus the 
fact that the finite aspect of space-time only appears in the evaluation of the number of 
allowed modes in the cavity by supposing that they are associated to TEM{y) waves with 
multiple perfect reflections along a given direction in the cavity. In fact, in order to correctly 
solve this problem, on must, along the chosen direction, use pairs of TEM waves of direction 
±9 n fulfilling limit condition (transversal resonance) with respect to the surface. This leads 
to developing the electromagnetic field in TE-TM modes, a set of complete and orthonormal 
functions for the quantization of a field in a box. 

By applying the theory of guided light propagation to the blackbody-resonator problem, 
we could introduce solutions in terms of finite field modes and also gain valuable insight into 
the dynamical processes allowing a blackbody to reach thermal equilibrium.The solution 
in finite spacetime with eigenvalues and eigenfunctions leads to the normalization of the 
exchange of energy between electrons in the volume S5 (S being the cavity area and 6 being 
the penetration depts of the EM field) and the EM modes oscillating in the total volume 
(S5 + V). The Planck constant ( h ) then appears as a normalization factor with respect 
to the characteristics of gold (as exemple). To the knowledge of the authors it is the first 
time an analysis of the Blackbody radiation including ’’power mode coupling” is presented. 
This mechanism, understandingly overlooked in Planck’s papers, also shine new light to 
the notion of decoherence in classical electrodynamics. The reader will find in table [I] a 
comparison between the approach proposed in this paper and the one generally presented 
in quantum physics. 
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Quantum Physics 



1. Plane wave development. 
(k 0 )(photon-like) 

l\k 0 \={kl+k;+kl) 112 ) 


2. Perfectly reflecting wall 
i.e. no loss, no coupling. 

3. Stationary wave at equilibrium. 

E{r,t)=q{t)E{r ) 

4. Genius empirical approach 
Energy=hco 

Planck-Einstein relation. 

5. Modal space- 
Number of modes 

, t/ \dv 8nv 2 dv 



6. Bose-Einstein statistics 
lead to Planck 
distribution. 


Modern electromagnetism 



1'. TE-TM mode development of 
EM field : complete set of 
orthonormal solutions. 


2'. Loss associated within volume 
(S5) 

3'. Coupling in the total 
volume (\/+S5) ^ Power 
Coupling Theory. 

4'.h=normalization factor from 
exchange energy (that paper) 


5'. In addition of 5, there 

is important coupling between 
high and low frequency 
modes as shown in that 
paper. 

6'. Planck's analysis use of 
Boltzmann statistics. 


TABLE I. Comparison between quantum physics analysis and modern EM analysis of Blackbody 
radiation. 

An analysis to be cited 0IB] based on zero-point radiation, pre-” suppose” the existence, at 
zero temperature, of’’classical homogeneous fluctuating radiation”, and as, now, in quantum 
physics with ’’quantum vacuum fluctuations”, use an ”ad hoc” correction to the definition 
of the equilibrium. 

In this article, we will first present a few reminders of relevant waveguiding topics: the 
expansion of an EM field into ideal fundamental modes, the EM fields in resonant monomode 
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metallic cavities with a modified analysis conserving a privileged propagation direction, and 
EM fields in multimode metallic cavities. This will introduce mode coupling as an essential 
feature of EM multimode guided propagation allowing fast power exchange between the 
modes. Second, we will turn to the study of blackbody resonators in the multimode regime 
where mode coupling, in addition to absorption and emission of radiation, is needed to reach 
thermal equilibrium. The average EM energy per mode in the cavity will be evaluated know¬ 
ing that when low-loss multimode guided propagation is in equilibrium, all modes possess the 
same attenuation coefficient. To relate the EM energy to the cavity temperature, optimal 
energy exchange between the modes and the cavity walls will be examined at maximum syn¬ 
chronization with the free electron mean relaxation frequency cu r , related to their mean free 
path between collisions. The reduced Planck constant (h) and the radiative energy density 
distribution will be shown to emerge naturally from this relation obtained from EM anal¬ 
ysis, h, becoming a derivable physical constant. The uncertainty principle is also emerging 
from this a priori classical theory. Furthermore, from a microscopic-macroscopic analysis, 
a frequency relation between the gold (reference) and others materials of a large group of 
good metallic conductors will be derived, thus establishing a link between the microscopic, 
cu T , matter oscillator (even magnetic one), and the macroscopic Planck distribution aspect 
of that problem. 


II. ELECTROMAGNETIC FIELD REPRESENTATION IN FUNDAMENTAL NOR¬ 
MAL MODES 

Guided, radiative and evanescent modes of waveguides constitute a complete and or¬ 
thonormal set [7]. Such a set can therefore form a basis for a series expansion of an EM 
field. Figure [l] represents the domains of these modes for a waveguide of refractive index 
ri\ in a medium of refractive index n 2 < rq, where /3 is the propagation constant of the 
observed mode.Without coupling, the power in each guided and radiative mode is trans¬ 
ported individually since /3 is real. In the case of evanescent modes having pure imaginary 
/3, a combination of progressive and retrograde evanescent modes is necessary to transport 
energy, they cannot do it individually um. 
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FIG. 1. Modal domains of a waveguide with refractive index ni in a surrounding medium with 
refractive index 77,2 < n\, where /3 r and are the real and imaginary part of the modal propagation 
constant respectively. 


One can expand the EM field of an irregular waveguide in terms of “ideal (normal) 
modes”. The terms of the expansion, then function of the propagation direction z, are so¬ 
lutions to a set of coupled differential equations. Another method exists to obtain a set of 
coupled differential equations describing mode coupling. One can formulate the EM field of 
the irregular waveguide at a point z a in terms of the modes of an hypothetical local waveg¬ 
uide. These modes are called “local modes”. However in this article, the first description 
based on ideal (normal) modes will be used as it is more appropriate for regular waveguides 
with ideal geometry and therefore lends itself well to the analysis of the blackbody problem. 
For a more detailed analysis of ideal mode coupling, the reader may refer to the third chapter 


of the book by D. Marcuse |H] which we will use in section III B of this paper. 


It is important to note that boundary conditions are not applied explicitly to obtain 
the solution that will be described. The ideal modes expansion method generates solutions 
that apply to all space, since the boundary conditions are contained in the mode equations 
composing the expansion. These boundary conditions are automatically verified by the 
solutions, thus there is no need to bring these relations later in the process. This modal 
expansion also speeds up numerical simulations in general by one or two orders of magnitude 
due to the finite number of modes comparatively to the spatiotemporal model which deals 
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with base equation of infinite dimensionality. 

III. RESONANT METALLIC CAVITIES 

A. Monomode metallic cavities in modern electromagnetism 

At first, we will follow an analysis that can be found in electromagnetism textbooks, 
see for example J. D. Jackson [iJj. and we will later adapt it to the type of analysis done 
in quantum physics. The EM field inside the metallic cavity will be expanded into ideal 
(normal) stationary modes, i.e. normal modes from short-circuited guided propagation 
along the z axis which allows us to maintain a preferential orientation. We start by choosing 
a metallic cavity of rectangular shape, in order to reproduce the case of ” quantization in a 
box”, with dimensions a, b and d, where d is the depth of the cavity along the 0 axis. 


X 


a 



Cutoff frequency for guided propagation along z: 



FIG. 2. Geometrical parameters of a rectangular metallic cavity and associated frequency relations. 
Since the z axis is designated as the axis of propagation and observation, = 2 d/n. 

For our modal expansion, we will label the mode in the x direction with m, in the y 
direction with p and in the z direction with n. Considering power losses in the cavity, one 
finds that the energy stored inside the cavity for a mode (■ mpn ) is given by the following 
equation: 



trans 


(i) 


7 












where 


F'trans 



7 r 2 c 2 

4u mp 


\^{x,y)\da 



( 2 ) 


with il?(x, ?/) representing the distribution of the transverse field associated to a ( mpn ) mode 
where n = 2d/\ g ^ n . 

At this step, authors often express the cavity volume V explicitly in order to determine 
the cavity’s Q-factor: 


Q = 




v 

— x (Geometrical Factor) 

DO 


(3) 


M conductor 

where V and S represent respectively the volume of the cavity and its surface area while 
S is the penetration depth of the field inside the walls of the cavity. The relationship (J2| is 
also used to determine the well-known density of modes 


N{y)dv/V = 8nv 2 dv/c 2 


(4) 


where N{y)dv represents the number of modes in the frequency interval [u, v + 5v\. Usually, 
one would obtain an energy density by solving an integral over only one frequency band 
which is also equivalent to using the cavity volume V explicitly. 

However, we will not make V explicit at the moment to keep the ’’behavior of time 
and space coupled”. We will continue our analysis with the short-circuited propagating 
wave approach, which maintains the analogy with the one-dimensional harmonic oscillator 
by preserving a preferential direction of propagation. Since d = nX gtTl /2, one can rewrite 
Eq. ([Tj) in the form: 

u A„ ,.u'/ . n/,i _ 

(5) 


j j _ ^^g,n^mpn jp _ mpn r, 

v-A n Z iA -l t.rn.n.fi t 


2,71 c 2 


trans 


y gr,n 


where v gr n is the group velocity of the guided short-circuited mode n in an empty cavity as 


defined by 


^ gr,n C\ I 1 


U, 


mp 


U, 


( 6 ) 


mpn 


and F trans is defined in Eq. (j2]). 

Important results are obtained by analyzing the field directionally along z for set m and 
p values. This step is essential to establish a link with quantum mechanics as the observable 
of the stored energy in the cavity can then be expressed as: 


TiUJrfJipjl 

U n = 

Vgr,n 


lx 


m P\TMTE 


-» 


nv, 


mpn 


lx 


m PiTM,TE 


(7) 










where the rightmost expression is derived for the asymptotic limit far from the cutoff fre¬ 
quency (uj mpn S> oj mp ) with Xm P designating a constant for TM and TE modes. 

If we assume lossless, perfectly reflective cavity walls, a def ini te link is thus established 
between modern electromagnetism and quantum physics through the quantization rw mpn 
of the energy. However, the effects of losses associated to finite conductivity are already 
included in uj mpn which is linked to the resonance frequency w 0 (mpn) of an ideal lossless cavity. 
Large values of ohmic losses are always causing a decrease in the resonance frequency, as 
seen in the following relation: 

cn-cno = Acu « (8) 

Thus our analysis does not make the typical, yet unrealistic, assumption of perfectly con¬ 
ducting walls. Let us add that when the system is left undisturbed, the energy initially 
stored in the mode n of a realistic cavity will decrease exponentially with a time constant r 
inversely proportional to Q, i.e. r oc [ojq/Q). The outward power flow is then equal to: 

Pn Vgr,nU n YUdmpn Ft.rans• (9) 

In the classical electromagnetic approach to the blackbody problem, the number of modes 
in the modal space was determined in conditions more restrictive than the asymptotic limit 
of the group velocity of all modes in the cavity tending toward c. A TEM wave method 
was used with v gr taken exactly equal to c, but the propagation of infinite TEM waves is 
forbidden mdividually inside a hollow metallic waveguide. How can one easily apply the 
finite boundary conditions of a blackbody in this situation while preserving the coupling 
between space and time variables? It is perhaps not surprising that these approximations 
led to the failure of classical electrodynamics for the blackbody problem. Interestingly in 
quantum physics, one implicitly makes the v gr = c assumption of classical electrodynamics 
when quantizing the energy according to E = nhv (Iw and the corresponding states | 'ifn)). 

In contrast, the blackbody boundary conditions are naturally satisfied by the TE, TM 
normal mode description from short-circuited guided propagation. A few observations read¬ 
ily follow from this description: 1) the bounded TE, TM modes link local phenomena to 
their environment inducing a related decoherence; 2) a wave packet formulation can be based 
on the group velocity v gr and the losses from the cavity environment; 3) the need for a renor¬ 
malization process is eliminated since the TE, TM modes inherently satisfy special relativity, 
and are a set of complete orthonormal solutions ; 4) this TE, TM modes picture is valid for 
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arbitrarily shaped cavities as long as their size is sufficiently greater than the wavelengths 
involved. Because of the latter, one can neglect a slow envelope modulation of the EM field 
in space as well as small perturbations near the cavity walls in the asymptotic limit v gr —>■ c. 
The spatially rapidly varying part of the EM field inside an arbitrarily shaped cavity then 
corresponds to the TE and TM modes of a rectangular cavity; each 3-dimensional period of 
this spatial variation will be henceforth referred to as TE-TM cell. In this case, all cavities 
have the same relation for the electromagnetic mode density which is also the one found in 
classical electromagnetism. 


B. Multimode guided propagation and mode coupling 

As was noted by Planck himself ra. a complete description of blackbody radiation re¬ 
quires another mechanism than light emission and absorption to explain how energy is redis¬ 
tributed between different frequencies to reach the thermal equilibrium distribution from an 
initial arbitrary distribution. In the present paper, we show that the mode coupling mech¬ 
anism introduced from guided EM theory in section [TT] has the necessary characteristics to 
ensure such an energy redistribution. We thus return to the EM field representation in ideal 
(normal) modes in the case of multimodal propagation with mode coupling. This involves an 
infinite set of integro-clifferential equations. However, this is an unsolvable problem except 
if certain conditions can be set to reduce its complexity. Marcuse demonstrated that it is 
possible to obtain a set of coupled equations bringing into play the average power in each 
mode instead of their amplitude |8]. This demonstration holds for waveguides of arbitrary 
length, as long as the coupling is weak. It is based on the experimental observation of the 
average power of a mode in highly multimodal propagation with weak coupling. After an 
adaptation length, this power averaged in z for a given mode v is of the form: 

pjjnoy) ( 2 ) = Aj/ exp(-az) (10) 

where cr is independent of v and will be found to represent the eigenvalues of a set of 
homogeneous equations. 

In the absence of coupling for multimodal guided propagation, an attenuation coefficient 
2 a u can be attributed to each uncoupled mode to represent the power flow carried by a 
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mode along the waveguide as follows: 


P u (z) = P u (0) exp(-2a u z) 


( 11 ) 


where P u ( 0) corresponds to the power flow of mode v at z = 0. The shape of the power 
distribution is stable only when the mode, or a small subset of modes, with smallest losses 
carries power while the other modes have disappeared completely. The total power carried 
by the set of all modes does not decrease exponentially in this case. 

In the presence of coupling induced by a mechanism independent of z, the power redis¬ 
tributes itself between the different modes, following the combined effects of individual losses 
and coupling exchanges. If this coupling is weak enough to keep losses very low, but strong 
enough to overcome individual losses, it tends to equalize the power carried by each mode 


and to establish a stable distribution of its power characterized by Eq. (10). The ratio of 
average power carried by each mode in relation to a reference mode then no longer depends 
on the length coordinate. The same attenuation coefficient can therefore be attributed to all 
modes.The weak mode coupling regime allowed Marcuse to obtain, using the perturbation 
theory, the following average power coupling relation: 

}(moy) N 

— p(moy) \( p(moy) _ j 


dPh 


dz 


= — 2 ( P.l m °“ ) - Pjr y) ) 


( 12 ) 


12=1 


where p and u indicate modes in this set and C^ represents the coupling coefficient as noted 


by Marcuse. If at a location z, only the mode p carries energy, (12) reduces to a simple and 
intuitive representation 


dPy 

dz 


N 


— | 2 I P^ 


12=1 


where the term of thermic loss is augmented (or compounded) by the sum of all the cou¬ 
pling coefficients indicating that the mode p loses power to all the other modes. In the 
complementary case of a mode p at location z not carrying any power, we then obtain from 


( 12 ) 


dP 

dz 


N 




12=1 


The positive derivative shows that the mode p gains in power at the detriment of all the other 


modes. A set of homogeneous equations with N unknown is obtained from (10) and (12). 
The N eigenvectors A^f 1 and associated eigenvalues result from the solutions to this set 
of equations he.. The evolution of the propagating held is represented in figure [3j 
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1 Adaptation zone 

Steady state distribution zone 
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FIG. 3. The initial power distribution at z = 0 enters an adaptation zone with a constructive chaotic 
power distribution up to the beginning of the steady state power distribution at z = L s t = 
where K is related to the tolerated signal-to-noise ratio. No physical meaning is associated to the 
eigenvalues a ^ and eigenvectors of order N superior to 1. 


The steady state power distribution obtained by Marcuse decreases exponentially with 
cd 1 ) being independent of v. 


Pjr oy \z) = 


N 




v= 1 


4 1} exp(-cr (1 ^) 


(13) 


for z > L st with the parameters represented in figure [3] L st defines the distance at which 
this steady state is reached with 


L s t — 


In k 


(14) 


(erf 2 ) — erf 1 )) 

where the value of k is adjusted according to experimental observations. The relative 
notion of stability (equilibrium) of information propagation expressed by the tolerance to 


a given level of noise is represented by relation (14) where n stands for the tolerated noise 
level. This relation shows that the equilibrium state is nondeterministic which is related to 
the fact that uncertainty relationships can be obtained from classical electromagnetism in a 


finite-spacetime as we will show in section V B After the adaptation zone 0 < z < L st , the 
distribution of power versus mode number depends only on the shape of the first eigenvector 
and is completely independent of the initial power distribution. Mode coupling thus 
tends to distribute the power evenly over all the modes as long as it does not introduce strong 
radiation losses. This condition is fully respected in the case of hollow metallic waveguides 
and also in low loss dielectric guides. 
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As can be expected for equilibrium systems, one cannot work its way backward from 
measures of the steady state to find out the initial power distribution over the modes (lack 
of reversibility). At best, Marcuse’s model allows us to know which initial state should 
be used to reduce or eliminate the adaptation zone as noted by Marcuse. 


IV. MULTIMODE RESONATOR AND THE BLACKBODY PROBLEM 


The typical approach to study the EM field in a closed cavity is to assume the system has 
reached its steady state in order to separate the time and space variables: E(r, t ) = q(t)E(f). 
However, by doing so we lose all the information about the transient state and modal 
coupling effects, since the relation between space and time is eliminated. The coupled 
power equations theory summarized above can address this issue and also ensures that the 
blackbody is an ergodic system even from a purely electromagnetic perspective. Indeed, 
the stationary states of very high order in a closed cavity are preserved as can be expected 
for a low-loss multimodal system (highly reflecting walls) with relatively important mode 
coupling redistributing the energy between the modes. These stationary states are mostly 
established by stable propagation inside the cavity while the radiation within the adaptation 
zone contributes noise to the system. 

Mode coupling was also confirmed experimentally m through the observation of fre¬ 


quency independent losses as predicted by the exponential power decay (13). For all modes 
circulating in a realistic cavity, we can therefore define the same characteristic length z i/ e 
traveled over ■// round-trips in a cavity of length d before they have lost 1/e of their power: 


Zi/e = (cr (1) ) 1 = r]d = 


^gr,n^~ri’ 


(15) 


The right hand side defines the corresponding lifetime r n of mode n in the cavity using the 
group velocity v gr>n from Eq. Q. We also note that the v gT)U T u product is independent of 
frequency and must be the same for all modes u. When u) mp c j mpn with the limit v gr ^ n — > c, 
the modal lifetime r n = t' also becomes frequency independent and the cavity’s Q-factor 
from Eq. (|3]) will be: 

Q = ujr '^—a)- ( 16 ) 

echo 

Independence of frequency thus extends to spectral broadening of the modes as well, since 
Au = u/Q ~ ccrd). Once equilibrium is reached, all the individual modes have the same 
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properties cr^v^Au and r', except for Q becoming a linear function of uj. 

To complete our analysis, the nature of the coupling between modes will be characterized. 
For the metallic blackbody cavity, ohmic losses are the dominating contribution to mode 
coupling. The Drude model tells us that absorption for each electromagnetic mode will be 
proportional to cu 1 / 2 in the microwave frequencies while in the far-infrared region, it tends to 
a constant for ca uj p , (see figure [4jb)the plasma frequency. In the near infrared however, 
the absorption function is significantly dependent on uj t , which is a damping frequency linked 
to the free electron relaxation time [T2] , For an isolated mode with resonance frequency uj 
very far from the cutoff frequency a i c — The modes with uj close to c u c do not obey this 
law, since their geometrical factor is different. Furthermore, the number of modes in the 
modal space no longer obeys the classical asymptotic relation. — losses in the cavity are 
also proportional to uj 1 ' 2 [9]. It is more convenient to represent the losses for uj uj c in 
function of uj, oj p , and uj t using the complex dielectric function e c as given by the Drude 
model US E]: 


e c — 6i + ie 2 — Coo 


0J„ 


uj 2 + iujuj r 


(17) 


where 


UJ p 

UJ T 


27rc\ 
1 

2ttct 


AnNe 2 

771 * 

cm -1 . 


cm 


(18) 

(19) 


In these equations, N is the free electron density, e is the charge of the electron and rn*, 
it’s effective mass. The high-frequency dielectric constant is represented by while the 
free electron relaxation lifetime is given by r. It should be noted that u T is evaluated from 
measurements of the complex dielectric function e c . 

Our analysis of the blackbody problem based on electromagnetism relies on the relation 
between uj t and the temperature. Finite held volumes in the cavity (electromagnetic en¬ 
ergy) and finite microscopic interaction volumes in the cavity walls (thermal energy) are 
exchanging energy by coupling. When the cavity sustaining the radiation held is kept at 
constant temperature, the best energy exchange will occur at the optimum synchronization 
frequency uj mpn = uj t . At this equilibrium, all modes travel the same number of round-trips 
and, in their asymptotic limit, have the same lifetime r’ in the cavity. We can then consider 
our resonant multimodal system with synchronization frequency uj t equivalent to a lossless 
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monomodal resonant system with the single resonant frequency u> 0 = ui T (see figure |4ja). For 
the mode at uj t and those nearby in frequency at a position z where p^ moy ^ is maximum, 


Eq. (12) yields: 


dP, 


(' moy) 


dz 


N 


E^- 


v=\ 


p(moy) 


( 20 ) 


at thermal equilibrium. The exchange of thermal to electromagnetic energy and vice versa 
is due to the coupling term T = C/w This energy exchange per unit of TE-TM cell 

in the propagation direction is expressed by equating Eq. ([ 7 ]) to the thermal energy as per 
the equipartition theorem: 

2 QfcT) = T ■ u; T = T • oj n . (21) 

The factor 2 comes from the transverse electric and magnetic fields contributing a degree of 
freedom each [T5j. As there are n cells in each independant propagation direction, the total 
energy in the mode of indice n is: 

E n = n uj n • T. (22) 

We can thus determine the value of the coupling constant T with a reliable value of lu t 
at normal temperature (293.15 K). Since the definition of a Kelvin jlbj (K) makes the 
temperature unit dependent on the property of a material, we must use a reference conductor 


satisfying the Drude model such as gold to ensure the accuracy of relation (21). Brandli and 


Sievers have shown that measurements of e\ and 62 for gold in the infrared follow the Drude 
model very well, without using adaptive parameters for the number of free electrons per atom 
as is often done for other materials mi- Other arguments in favor of using the oj t value 
of gold are its very good experimental agreement for the conductivity according to Ohm’s 
law and negligible interband contributions to its absorption in the infrared. The average 
value measured for gold with different methods is (u T ) gold = 2.11 x 10 2 cm" - 1 [131 HH [Ill- 


According to ( 21 ), we then obtain the constant value: 


Ty 0 id = 1.02 x 10“ M J ■ s = h. 


(23) 


This result agrees well with the quantum mechanical analysis. Here, the Planck constant 
appears as a normalization factor between electromagnetic and thermal energy exchanges 
at thermal equilibrium. Reuniting these two domains, this result expands the signification 
of the Planck constant beyond the commutator function it fulfills in the quantization of a 
particle-electromagnetic field system. Electromagnetic coupling therefore plays an essential 
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role: it enables the cavity to reach the absolute stable equilibrium required by the second 
law of thermodynamics. This equilibrium is achieved by minimizing cavity losses with 
power transfer from modes having high losses to modes having lower losses as well as with 
the well-known emission and absorption of radiation. To resume the ” stationnary” states, 
equilibrium, is mostly established by stable propagation inside the cavity while the radiation 
within the adaptation zone contributes to noise of the system, as represented by the short- 
circuited propagation model of figure [3j 


a) 




FIG. 4. Schematic distribution for a) a monomodal system with u 0 > ui m ax , b) losses of isolated 

modes (—) and loss cd 1 ) at equilibrium (-) for a multimodal system with representation of 

the two initial groups of modes in term of their resistance, c) resultant energy distribution at 
equilibrium where u T plays the role of uj 0 and u T < co ma x ■ 
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Two different situations can occur in the process (figure [4]). In the first case (4.a) where 
the losses are important, in a closed or quasi-closed system, it causes strong dissipative 
power transfert from high frequency modes (of higher loss) to lower frequency modes (of 
lower loss), and ultimately, it brings the collapse, at the equilibrium, in a few modes or 
only in one mode as in figure |4ja), where u Q corresponds to the resonnant frequency of an 
equivalent system without loss, and where uj 0 > oj max . In the second case, for a multimodal 
system with very low loss and relatively strong coupling, as, for example, in metallic thermo¬ 
static cavities, the high frequency modes (higher loss), mainly, release power by intermodal 
coupling to low frequency modes (lower loss) and u T , the relaxation frequency, of matter 
oscillators, plays an equivalent role as u Q , but now with uj t < u> m ax ■ Mode coupling thus 
tends to distribute the power evenly over all modes as long as it does not introduce strong 
radiation loss. That condition being respected in the cases of hollowo metallic resonator and 
in low loss dielectric resonators. It should be noted that in the first stage of the equilibrium 
process, modes of two frequency ’’zones” contribute to the coupling in different manners 
(see figure 0b); in zone B (u n > ui T ) the resistivity of an equivalent system is initially 
positive, while in zone A (oj n < U) T ) the resistivity is initially negative, meaning that its 
contributions modified the held without important contribution to p op (optical resistivity). 
In a way, the modes act, in part, as in a semi-conductor material and thus, our analysis 
presents an intuitive approach to the study of semi-conductor optical cavities following the 
finite-spacetime paradigm in classical electrodynamics. That research is out of the scope of 
that paper. 

Also of importance, is the following remark about the ’’dissimilarity” of the modes related 
to finite-spacetime through their oj c , cutoff frequency, function of the mode order, and their 
distinct coupling effects. So it is not appropriate, in a finite multimodal system, to assert 
that ’’all elements are alike” and that ’’all dissipative effects” act to lower the maximum 
frequency of the resulting energy density spectral distribution. Modes of frequency u> n < co T 
act to enhance, by multimodal coupling, the uj max of the spectral energy density distribution 
(Planck distribution). That effect is essential to the establishment of the equilibrium in the 
system. These two characteristics are essential to explain how the equilibrium is attainable 
within the finite-spacetime paradigm and the TE-TM held developpment based on a complete 
set of orthonormal functions. Otherwise, as in the space-time paradigm where the photons 
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are all alike, one must accept an ”ad hoc” commutator function for h, uncertainty relations as 
principles and ” quantum vacuum fluctuations” as derived from h and uncertainty relations, 
that notion being used to explain the Casimir effect [ID]. 


A. Planck’s law 

The previous results obtained from a classical electromagnetic analysis demonstrate that 
the equilibrium state in the cavity is reached and maintained through hv quanta character¬ 
izing exchanges of thermodynamic and electrodynamic energy. Each exchange of a quanta 
can be interpreted as a temporally finite excitation of a mode and in blackbody radiation 
conditions, the total number of excitations easily falls into the limit of large numbers. Since 
an electron can produce over time an arbitrary number of excitations fluctuating in energy, 
hence possibly but not always coupling to the same mode, the time-averaged energy dis¬ 
tribution of these excitations will follow the distribution of energy density as derived by 
Planck [2D], as example see page 19 of J2TJ. With the average energy of this distribution 
and the well-known density of modes N(u)du/V = Rn^dv / c 3 , which can be calculated from 
either the modal approach followed here or the classical counting of standing waves in the 
cavity, we express Planck’s law as follows: 

8n hv 3 1 huj 3 1 

p(u)du ——^- 7 — 77 =t- dv o- ——— 7 —77777 - du = p(uj)doj. 24 

C 3 gold TT^C 3 (X'gold ^ 

We chose to express the familiar relation for radiative energy density p(v)dv with the damp¬ 
ing frequency c o r (T) instead of the usual h/kT factor to emphasize the connection between 
temperature and materials. In the next section, we will analyze the relevance of other ma¬ 
terials, than gold, by obtaining a parametric frequential relation linking all good conductors 
to Planck distribution. 


B. Good conductors metallic cavities 


A comparison against the gold must be established to obtain a radiative energy density 
distribution p(co)dco obeying Planck’s law as observed for several metals. This is done by 


introducing a comparison factor Gx in the maximum synchronization relation ( 21 ) for an 
arbitrary metal labelled X, i.e., 
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GxkT = Constant ■ (c u mp n)x = Constant ■ ( u T )x ■ (25) 

The factor Gx allows the scaling of temperature in other materials against the gold 
(as example) with kT defined as h(u> T ) go i d we obtain Gx = (w T )x/(w T )goid from the above 
relation. In the time domain, we can interpret this ratio as the relative number of cycles per 
unit of time over which an energy quanta is exchanged, which is consistent with the Planck 
constant units in joules per hertz. In the spatial domain, this corresponds to the ratio of 
TE-TM cells per unit of length over which coupling takes place, in relation with the weak 
losses and penetration depth of the electromagnetic field in the cavity walls. Additionally, 
different material will have a different density of electrons emitting and absorbing radiation, 
which will be included in Gx via the ratio of free electron concentrations N e . We must also 
include materials parameters influencing the overall system responds through the synchro¬ 
nization between the electromagnetic modal electronic and the damping frequencies. 


The microscopic model of the Ohm’s law implies (18J that the courant-voltage relation 
relation at the input-output of a material is due to the fact that the associated electric field 
adds a small (regular) speed to the chaotic one of the free electrons in the material. Taking 
into account relation (19), one finds that [TB] the relation frequency is 


where 


w T (cm 1 ) 


N e e 2 p 

m e 27rc 


(26) 


N e = electron number/cm 3 
e = electric charge 


p = resistivity 


rn P = electron mass. 


That relation (26) applies only to materials corresponding to the Drude model, as for 


example, gold and silver. At normal temperature (293 K) their calculated results from (26) 
are (co t )a u = 2.11 x 10 2 cm _1 and {uj tau )Ag = 1-42 x 10 2 cm“ 1 while their mesured ones are 
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(u t )au = 2.15 x 10 2 cm 1 and (a k)aq = 1-45 x 10 2 cm 1 from [IT3l \li\ . 


So, in the case of these materials, the parametric frequential relation Gx is reduced to 
only ratios of electronic density and resistances; i.e. 


(w T )x _ (N e ) x px 

(^t) gold (A^ e ) gold Pgold, 


(27) 


The number of the materials obeying to relation (J2T]) is quite small, Al already departing 
of it by almost 13%. The authors of reference [T3( [T3j note, without surprise, that at high 
frequencies, Au, Ag and Al, follow well the free electron model, but that the overall pic¬ 
ture does change for the others materials (interband and surface effects). Furthermore, in 
a metallic thermostatic cavities, at the microscopic level, the resulting resistivity is mainly 
determined by the modal high frequency phenomena, those modes releasing most of their 
power by coupling to the lower frequency modes (see figure [4]). The dominant physical mech¬ 
anism in our approach is ultra fast coupling compared to slower absorption and re-emission. 
Hence, one have to choose appropriate effective N e ,m e ,p. In classical electrodynamics, a 
’’good conductor can be modelled as a perfect conductor with the idealized surface courant 
replaced by an equivalent surface courant, which is actually distributed throughout a very 
small, but finite thickness at the surface. Thus, the surface resistance has to replace the 
volume resistance; i.e. 1/a —* 1 jab, where a is the conductivity, and S the skin depth” [9j 
chap.8]. In a preceding section [Tv] it was shown that the ’’modal” quality factor, in multi¬ 
modal cavity resonator, was, 


( Q)x - 


UJn 


ca 


(i)’ 

x 


(28) 


where, here, the subscript X is for a chosen material. Remembering that, at equilibrium, 
’’a*- 1 )” and ’’r'” are frequency independant function, then (26) can be written 


cu T (cm *) 


m* e 6 2i tc 


(29) 


where p' is in H • cm, p' % p, N e , the electron density (free and valence electrons) and m* e 
the effective masses of these electrons in their inertial opposition to the movement. 

As from reference [9j, cited here as equation ([3]) 
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Q = 


p 


— ) x (Geometrical Factor) 


pconductor \ S8 / 

Using equations (28) and ([3]), and the plasma frequency of the material X; i.e. 


@i 


C u p )x - 


27rc 


AirNpe 2 


ml 


1/2 


(30) 


we can, then, obtain, in substituting (30) in (29), and from (J3]) and (28) 


a {1) - 
° x ~ 


p'xH 


IX 


S 


(p 


c)X 


2 V ■ ( geo.fac.) /j ^ ^ 

As following [9], the ’’geometrical factor” evolutes towards a constant for hight mode 
number, the term S /{V -(geo.fac.)) appears as a constant for a cavity (resonator) of specified 


form in the analysis for different materials. Relative to gold, equation (31) implies 


G x = 


K) 


a 


(i) 

x 


_ _ (Vp)x (p c )x p'x 

t)au &Au ( W p)Au (pc) Au Pau 


(32) 


Thus, the a ^ must be measured in an equivalent metallic guide in equilibrium propaga¬ 


tion and the p' x obtained from (32). The phenomena of multimodal coupling being related 
to a strong power exchange from high frequency modes to lower frequency modes, the p' x 
will, then, correspond to the p op ticai of references [EiEj, in the case where (jf^j — 1- In 
other cases, as Fe, Ni,Co, the result has to be labelled (p' x , p Cx )optical- As for Au and Ag, 
their p op u C ai - Po, (u p )au - (u p )a 9 and their (j^j ~ 1, @ becomes 


r (1) 

r Ag 


°aI PAu 


[U. 


~) Ag 


(^t) Au 


(33) 


Table [TT] shows the results obtained for the normalization of an ensemble of good conduc¬ 
tors as based on references [T3] and [IT] , from where the different measured values of (ca T ) x , 
gold assuming a role of reference at a T = 293.75K, i.e. (oj t )au = 2.04cm _1 . It shows also 
that relative coefficients G x are in ’’accordance” to the ratio (uj t )xI(u t )au and of cr^/cr^. 
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That analysis shows that a good conductor metallic cavity, when maintained at constant 
temperature, acts, so far as its spectral density distribution is concerned, as if the field 
as equilibrium is only ” confined” within the volume ”V”, and h being a physical constant 
related to a normalization in its development in a complete set of orthonormal functions. 
Thus the overlooked finite-spacetime paradigm in blaekbody radiation analysis leads to 
establish the ” ad hoc” commutor function of h in quantum physics, in this case, on a sound 
physical theory, the Maxwell theory. Moreover, it explains the statistical aspect of the 
vacuum field through multimodal coupling of two kinds, while in the restricted spacetime 
paradigm of the quantum physics, that aspect rest upon the notion of ’’quantum vacuum 
fluctuation”. Thus, one cannot assert that the notion of decoherence by multimodal coupling 
in classical electrodynamics is emerging, in this case, from quantum physics. It is the result 
of finite-spacetime analysis of electromagnetic fields in multimodal modes and states and 
the ’’smearing off’ of informations that can be obtained from the systems. 


V. CONCLUSION AND OUTLOOK 

Connections between quantum mechanics and electromagnetism have been established 
throughout the analysis of the blaekbody problem presented in this article. We came to the 
conclusion that coupling between modes of the cavity, following the coupled power theory 
for highly multimodal guided propagation, is an essential element in the analysis of the 
blaekbody cavity and sheds new light on the physical interpretation of this century-old 
problem. The main features of the approach followed were: i) replacing the basis of infinite 
plane waves by finite TE, TM modes to keep the system of electromagnetic differential 
equations coupled with space and time, ii) analyzing the interaction of the EM field with the 
cavity walls based on the temporal lattice of damping frequencies created by the free electron 
relaxation lifetime, and iii) considering the maximum energy exchange in this interaction 
occurring at synchronization between electronic damping frequencies and electromagnetic 
modal frequencies. The Planck constant naturally appears at the outcome of this analysis 
using material parameters from gold as example, the ’’path towards equilibrium” providing 
a classical value for h. We have also shown that the same process applies to a large group of 
good conductors metallic cavities. It follows that the Planck constant can now be derived as 
a normalization factor between thermal and electromagnetic energy exchanges at thermal 
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equilibrium and the ’’ultraviolet catastrophe” is no longer an unexplained phenomenon in 
classical electrodynamics and it shines new light on the notion of decoherence. We provide 
here only a brief outlook on some other basic issues. 


A. Fluctuations 

Our focus here was on the mechanism proposed to reach equilibrium in the cavity and 
the Boltzmann distribution has been the only result used from statistical physics. We how¬ 
ever want to emphasize that statistical mechanics is still fully applicable to the blackbody 
cavity problem governed by coupled power theory. This theory encompasses power fluc¬ 
tuations for an electromagnetic mode, or a subset of a few modes, away from the average 
power distribution of the entire set, hence preserving the stochastic nature of the system. 
Additionally, mode coupling enables a fast redistribution of power between them, hence ac¬ 
celerating the evolution of the system toward equilibrium. As stated previously, if the initial 
power distribution is unknown, it is not possible to apply a reversibility principle to such 
systems. 


B. Uncertainty principle 


Observable quantities of TE-TM waves as derived here from the coupled power theory 
and experimental measurements will obey the general recognized uncertainty properties for 
all EM waves: 


AvAt A — , 
47r 

AxAk A — . 
47r 


(34) 

(35) 


Relation (34) is readily seen to correspond to Heisenberg’s uncertainty principle AEAt ^ 
h/2 since we also derived E = hu as the quantization of thermal and electromagnetic energy 
exchanges. This suggests to interpret the uncertainty principle as a consequence of the 
finite space-time for the EM waves: without the infinite plane wave, there is no more pure 
frequency or wavelength, each finite wave excitation is associated to a frequency bandwidth 
and finite extent in space. It is possible to go beyond this first order approximation by taking 
into account the group velocity v gr of TE-TM waves in the the AxAk inequality through 
the mode momentum p = E jv gr = hv/cyjl — w^/w 2 where c o c is the cutoff frequency. 


24 




C. Photoelectric effect 


The photoelectric effect has traditionally been one of main justifications for the photon 
hypothesis to represent electromagnetic radiation. However in more recent years, it was 
shown that a semi-classical theory where only the photoelectric detector is quantized in terms 
of h embodies all the features of the photoelectric effect m In this article, we presented 
further evidence that the appearance of h in light absorption, emission and coupling does 
not call for the introduction of the photon concept. Also, temperature dependence of the 
photoelectric threshold is more directly linked to material properties in our analysis by 
using the damping frequency. Finally, we emphasize that the dominant physical mechanism 
in our approach is ultra fast coupling compared to slower absorption and re-emission, hence 
matching characteristic time scales of the photoelectric effect. 


VI. ACKNOWLEDGMENTS 

This work was supported by the Natural Sciences and Engineering Research Council of 
Canada (NSERC). We also acknowledge the precious help of Gabriellc Theriault, Alexan¬ 
dre April, Jean-Raphael Carrier, Pierre Tremblay and the recent collaboration of Anthony 
Richard. 


[1] G. Kirchhoff, The London, Edinburgh and Dublin Philosophical Magazine and Journal of 
Science 4, 20, 1 (1860). 

[2] W. E. Lamb, Applied Physics B 60 , 77 (1995). 

[3] G. Beck, Zeitschrift frir Physik 41 , 443 (1927). 

[4] M. Y. Han, A Story of Light, edited by W. S. P. C. P. Ltd (2004). 

[5] T. H. Boyer, Phys. Rev. 182 , 1374 (1969). 

[6] T. H. Boyer, Foundations of Physics 42 , 595 (2012). 

[7] H. Kogelnik, Theory of Dielectric Waveguides in Integrated Optics, edited by T. Tarnis 
(Springer-Verlag, New York, 1979). 

[8] D. Marcuse, Theory of Dielectric Optical Waveguides, 2nd ed. (Academic Press, Boston, 1991). 

[9] J. D. Jackson, Classical Electrodynamics, 3rd ed. (John Wiley Sz Sons, Inc., New York, 1998). 


25 



[10] T. S. Kuhn, Blackbody Theory and the Quantum Discontinuity, 1894-1912, 2nd ed. (University 
of Chicago Press, Chicago, 1987). 

[11] D. Marcuse, Principles of Optical Fiber Measurements (Academic Press, New York, 1981). 

[12] H. J. W. Miiller-Kirsten, Electrodynamics: An Introduction Including Quantum Effects, 1st 
ed. (World Scientific Publishing Co., New Jersey, 2004). 

[13] M. A. Ordal, L. L. Long, R. J. Bell, S. E. Bell, R. R. Bell, R. W. Alexander, Jr., and C. A. 
Ward, Applied Optics 22 , 1099 (1983). 

[14] M. A. Ordal, R. J. Bell, R. W. Alexander, Jr., L. L. Long, and M. R. Querry, Applied Optics 
24 , 4493 (1985). 

[15] W. H. Louisell, Quantum Statistical Properties of Radiation (John Wiley InterScience, New 
York, 1973). 

[16] H. Preston-Thomas, Metrologia 27 , 3 (1990). 

[17] G. Brandli and A. J. Sievers, Physical Review B 5, 3550 (1972), 

[18] N. W. Ashcroft and N. D. Mermin, Solid State Physics, 1st ed. (Brooks/Cole, 1976). 

[19] R. L. Jaffe, Phys. Rev. D 72 , 021301 (2005), 

[20] M. Planck, Eight Lectures on Theoretical Physics, Dover Books on Physics (Dover Publica¬ 
tions, 2012). 

[21] R. Eisberg and R. Resnick, Quantum physics of atoms, molecules, solids, nuclei, and particles 
(Wiley, 1974). 

[22] R. Loudon, The Quantum Theory of Light, 3rd ed. (Oxford University Press, 2000). 


26 








